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Abstract We report computer simulations of the interaction of seeding parti¬ 
cles with quantized vortices and with the normal fluid flow in thermal counter¬ 
flow of superfluid 4 He. We show that if the particle concentration is too large, 
the vortex tangle is significantly affected, posing problems in the interpreta¬ 
tion of visualization experiments. The main effects are an increase in vortex 
line density and a change in polarization of the vortex tangle, caused by the 
action of the Stokes drag of the viscous normal fluid on the trapped parti¬ 
cles. We argue that in the case of large particle concentration, typically used 
for the particle image velocimetry technique, the tangle properties might be¬ 
come significantly changed. On the contrary, the particle tracking velocimetry 
technique that uses smaller particle concentration should not be appreciably 
affected. 


1 Introduction 

Flow visualization is one of the most valuable experimental tools in fluid dy¬ 
namics. A variety of techniques exists which seed the fluid with contrasting 
agent that can be easily observed, e.g., by a sensitive camera [T]. Visualization 
techniques have already been used at low temperatures to study both classical 
and quantum cryogenic flows, even though their application is difficult, for 
various technical reasons, such as optical access to the experimental volume or 
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choice of suitable seeding particles. As for quantum flows and turbulence [2J (so 
far of superfluid 4 He, known as He II), application of visualization techniques 
face additional fundamental difficulties, due to the existence of two velocity 
fields and the interaction of seeding particles with quantized vortices [3]. De¬ 
spite these problems, the implementation of visualization methods in He II has 
led to the direct visualization of quantized vortices in |I] and important re¬ 
sults [5] about Kelvin waves and vortex reconnections [7], vortical structures 
around a cylinder 18), non-Gaussian velocity statistics [911101 and non-classical 
acceleration statistics El , crossover between quantum and classical behaviour 
[TO; particles’ trapping mechanisms m , and added mass effects El- In 
comparison with classical fluid dynamics, visualization of quantum flows is not 
yet firmly established and the interpretation of experiments poses important 
fundamental questions. 

The questions arise because He-II differs from classical liquids in several im¬ 
portant aspects [2]. For 1 K < T <T\, where visualization methods are usually 
applied, He-II can be described as consisting of two fluids - the inviscid super¬ 
fluid component carrying no entropy, and the viscous normal fluid behaving 
approximately as an ordinary Navier-Stokes fluid. Turbulence in the superfluid 
component can only exist in the form of a complex tangle of quantized vortices 
thin topological defects around which the circulation is restricted to single 
quantum of circulation k = h/M where h is the Planck constant and M is the 
mass of 4 He atom. 

This complex nature of He-II poses challenges for the interpretation of vi¬ 
sualization experiments. Here we consider potential problems that might occur 
in the interpretation of particle tracking velocimetry (PTV) and/or particle 
image velocimetry (PIV) techniques, both relying on observing seeding par¬ 
ticles (in recent experiments mostly solid hydrogen or deuterium flakes 0) 
suspended in the flow. These particles interact with both the normal and su¬ 
perfluid velocity fields [31I51HT] and can become trapped on the cores of quan¬ 
tized vortices. However, in most cases they are treated as non-intrusive, passive 
probes of the flow under study. It is therefore of great interest to try to assess 
the degree of non-ideality of the particles, which might lead to a distorted 
physical information about the quantum flow under study and, consequently, 
to a misleading conclusion on some important aspects of quantum turbulence. 

To this end, we perform a series of numerical experiments, extending the 
work of Mineda et al. M, which simulates the movement of seeding particles 
in the velocity field due to the counterflow tangle of quantized vortices. We find 
that a trapped particle deforms the vortex on which it is trapped, stretching it 
in the direction of the flowing normal fluid, via the action of the Stokes drag. 
A large number of particles can increase the vortex line density by up to 100% 
and change its polarization. We compare these more realistic results with those 
obtained by modelling trapped particles as ideal tracers of the vortices and 
find significant differences, both in the probability of trapping and the velocity 
statistics. 
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2 Computational setup and results 

We perform a vortex filament simulation of counterflowing He-II in periodic 
boundary conditions. The normal fluid velocity v n is uniform and statically 
prescribed. To a statistically converged counterflow tangle, we add inertial 
particles at random positions and initially zero velocities. The movement of 
both vortices and particles is described below. 

The quantized vortices in He-II are modelled, following the seminal 
work of Schwarz m , as one-dimensional spatial curves of arbitrary shape. 
These vortex lines (labelled s(£), where £ is the arc length along the line) 
induce a superfluid velocity v s given by the standard mmm Biot-Savart 
integral 



(1) 


with C denoting the entire configuration of lines in the vortex tangle. The 
movement of the vortices themselves is determined from the balance of forces, 
namely Magnus force and mutual friction, under the assumption that the 
vortices are massless (e.g., M)- The resulting equation of motion is 


s = v' s (s) + as' x (v n - v s ) + a's' x [s' x (v n - u s )], (2) 


where the dot denotes the time derivative and the prime denotes derivative 
with respect to the arc-length. The prime on the superfluid velocity in ED 
denotes the standard de-singularization of the Biot-Savart integral 

Particle dynamics. Several past studies were concerned with modelling 
inertial particles that are either far away the from cores of the quantized vor¬ 
tices dll or remain trapped on them m- Studies that address the full range 
of possible interactions of finite spherical particles and quantized vortex lines 
have already been performed [19], however, the computational complexity of 
the methods used there prevents the scaling of the simulation to high densities 
of vortex tangle and/or high number of particles. 

In the present study we adopt a mixed approach. Unlike EH , particles are 
considered as point-like objects. Particles sufficiently far from the vortex lines 
are considered free and interact with vortices only through inertial forces - 
following |18j , such free particles interact with both the normal and superfluid 
component inertially and with normal component also viscously so that the 
equation of motion for the the free particles’ velocity v p is 


d v p _ p s D v s As Du, 
df p Dt p Dt 


(3) 


T 


where Dv/D t = dv/dt + (v ■ V)v and r = 2a 2 p/9p, n is the viscous relaxation 
time, with p being the density of fluid or the particles (particles are assumed 
to be neutrally buoyant), a is the radius of the particles and p n denotes the 
dynamic viscosity of the normal fluid component. Time-independent and uni¬ 
form normal fluid velocity is used in the present simulations, therefore the 
second term on the right-hand side of m is identically zero. 
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Upon sufficiently close approach of a particle to a vortex, the particle can 
become trapped. Trapped particle experiences additional forces acting on it - 
namely the vortex tension, the Magnus force and the mutual friction |14l . It 
should be noted, however, that the particles are still assumed to be infinitesi¬ 
mally small and therefore a particle trapped on a vortex does not modify the 
superfluid velocity induced by the vortex, except possibly deforming the vortex 
itself. The modified equation of motion for the trapped particles becomes 

dUp _ v n - v p Ps Vv's 
df t p D t 

+ - s'_ ) + p s ns' x ( v p - v s ) (4) 

+ ( 70 s' x [s' x (v p -«„)]+ 7o s ' x ( v p ~ Vn )) > 

where T 0 = p s k 2 /Att \og[2y/l + l-/y/et;o] is the vortex tension (energy of the 
vortex per unit length) with l + , Z_ being the distances to the neighbouring dis¬ 
cretisation points along the line and £o ~ lA being the vortex core parameter, 
7 o and 7 q are mutual friction parameters related to the more standard a and 
a! [ 20 ], s'_ and s' + are right and left tangents of the vortex line (the line could 
be non-smooth at the site where the particle is trapped); s' is the tangent at 
the trapping site calculated as if the line were smooth, M e // = 3/2 pV p is the 
effective mass of the particle (V p being its volume) and At; is the maximum 
discretisation distance along the vortex, I. 6 xl 0 -5 m. 

From the point of view of the vortex, the equation of motion is changed 
from 0 for the single discretisation point that hosts the trapped particle, to 
Even though the equation of motion for only a single point is changed, 
this affects the local geometry of the line thus affecting the neighbourhood of 
this point as well. 

The present work differs from that of Mineda et. al m in that we consider 
both free and trapped particles simultaneously. This requires that the particles, 
once trapped, can de-trap and become free again. Algorithmically, the particles 
de-trap only if they are attached to a small vortex loop that is annihilated as 
a part of small-scale numerical cutoff. This mechanism roughly corresponds to 
a physical scenario (see Fig. [l]i) where a vortex loop collapses into the surface 
of the particle. An additional scenario of de-trapping is also captured (see 
Fig. Eb) - the trapped particle produces a cusp on the line sharp enough to 
cause reconnection producing a small loop that is immediately annihilated. 

This treatment of trapped particles is compared with the limiting case of 
particle so light that they do not affect the movement of the vortices in any 
way. In this case the particle is simply attached to a discretisation point along 
the vortex line and follows its movement - that is, its movement is described 
by ©. For the lack of a better name, this case will be henceforth referred to 
as “ideal” particle and the case described above as “non-ideal”. 

We performed several simulations with various numbers of particles N p 
(100, 250, 500, 750, and 1000) in a 1 mm 3 cubic computational box with peri¬ 
odic boundary conditions. All particles experience Stokes drag corresponding 
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Fig. 1 Two mechanisms of de-trapping 


to a sphere with 5 /tm radius. The imposed flow condition was counterflow 
with uniform, stationary normal fluid velocity v n = 0.55 cm/s and stationary 
imposed superfluid velocity. The initial condition of each simulation was a vor¬ 
tex tangle converged to a statistically stationary state with total vortex line 
length of about 6.7 cm. The initial condition for the particles was zero veloc¬ 
ity and random position, all the particles being free. All necessary parameters 
were taken from 23] at T = 1.9 K. 

After introducing the particles into the tangle, we observe two major 
changes in its properties. Firstly, as can be seen in Fig. [3 the particles in¬ 
crease the tangle density and for the highest number of particles this increase 
is nearly 100%. Secondly, the polarization of the tangle is also affected. Coun¬ 
terflow tangles are always partially polarized in the sense that the vortex line 
length in the transverse direction is larger than in the stream-wise direction, 
due to expansion of the favourably oriented vortex loops by mutual friction. 
Introduction of the particles is seen, in Fig. [3 to decrease this anisotropy - 
that is, to increase the relative length of the vortex tangle projected onto the 
stream-wise direction. 

Statistical properties of the particles’ motion, in the form of probability 
density functions (PDFs) of the instantaneous velocity component of free and 
trapped particles, are shown in Fig.[3for both non-ideal and ideal particles. In 
Fig. |3 the percentage of the trapped particles is also shown. Neither free nor 
trapped particles significantly dominate and, indeed, the PDFs feature con¬ 
tributions from both. The narrow peaks correspond to the uniform motion of 
the free particles determined mostly by the Stokes drag. The wide background 
distribution is a result of more chaotic movement of the trapped particles. 


3 Discussion and Conclusions 

Flow visualization using frozen hydrogen/deuterium particles is an important 
experimental tool of contemporary quantum turbulence research. Understand¬ 
ing the detailed behaviour of the particles is therefore of utmost importance for 
extracting information from the particles’ observed movement. Since detailed 
numerical simulation such as Ref. m is at this point prohibitively compu¬ 
tationally expensive for large number of particles and/or dense tangles, our 
attempt at providing some insight into particle dynamics is through the sim- 
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time (s) time (s) 

Fig. 2 (Left) Total vortex line length in 1 mm 3 periodic computational box for different 
particle number densities. (Right) Polarization of the vortex tangle, shown here as a ratio 
of the vortex length projected in stream-wise direction (a;) to that in transverse direction 
(mean of y and z). 



Fig. 3 Probability density functions of particles’ velocity component, for both non-ideal 
(right) and ideal (left) particles. Stream-wise (first row) and transverse (second row) com¬ 
ponents are shown. 
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Fig. 4 The percentage of trapped particles. Notice that the vortex stretching caused by the 
particles decreases the number of trapped particles by nearly a half. Notice also that the final 
ratio of trapped particles is nearly independent of the total number of particles, meaning 
that global changes to the structure of the tangle do not significantly affect de-trapping. 


plified point-like model used by Mineda et al. Mi but including trapping and 
de-trapping of the particles on the cores of quantized vortices. 

Our data show that the intricate interaction of the particles with quan¬ 
tized vortices can indeed become significant, and particles can cease to act as 
passive probes. This is evidenced by a significant increase in the vortex line 
density and by a change of the polarization of the tangle (Fig. [2]). These global 
tangle properties, however, are significantly affected for high enough particle 
concentration only. Careful particle tracking velocimetry (PTV) experiments 
that use low particle concentrations (typically about 1 mm~ 3 ) should not be 
appreciably affected. Results of methods that use higher particle concentra¬ 
tion, such as particle image velocimetry (PIV), should be interpreted with 
care. 

Moreover, one should remain cautious even if the low particle concentration 
does not affect global mean properties of the tangle. From the trapping rate 
(Fig. 01 one can see that both trapped and free particles give significant con¬ 
tributions to the overall statistical properties of the particle motion. Trapped 
particles, however, do not sample the vortex as it would have moved without 
this trapped particle - it locally deforms it and causes it to move with slightly 
artificial velocity. 

The effect of the back-interaction of the particles on the vortices on the 
statistics of the particles’ velocities can be seen in Fig. [3] It should be noted 
that the sharp peak corresponding to the free particles in the PDFs is generally 
not observed in the experiments. The most probable causes for this discrep¬ 
ancy, besides the generally wide distribution of particles sizes and shapes, is 
that the normal fluid component in the experiments is turbulent, while in the 
simulations it is laminar. 

The program used in the simulations is a modified version of the qvort 
code created by Andrew Baggaley and others, to whom we are grateful. The 
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